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We present results for the branching fractions and charge asymmetries in B ± -> h ± n° (where h 
— 7T , K ) and a search for the decay B° — ► 7r (, 7r° using a sample of approximately 88 million BB 
pairs collected by the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. We 
measure B(B ± — > ir n° ) = (5.5lo_g ± 0.6) x 10~ 6 , where the first error is statistical and the second 
is systematic. The B — > ^Tr" signal has a significance of 7.7a including systematic uncertainties. 
We simultaneously measure the A' ± 7r branching fraction to be B(B ± — ► K ± n°) = (12. 81^ ±1.0) x 
10~ 6 . The charge asymmetries are A„± n o = -0.031°;^ ± 0.02 and A K ± 7r o = -0.09 ± 0.09 ± 0.01. 



We place a 90% confidence-level upper limit on the branching fraction B(B° — > 7r°7r°) of 3.6 x 10" 
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The study of B meson decays into charmless hadronic 
final states plays an important role in the understand- 
ing of CP violation in the B system. In the Standard 
Model, CP violation arises from a single complex phase 
in the Cabibbo-Kobayashi-Maskawa quark-mixing ma- 
trix Vij 0|. Measurements of the time-dependent CP- 
violating asymmetry in the B° — > ir + ir~ decay mode by 
the BABAR and Belle Q collaborations provide in- 
formation on the angle a = arg [— ^d^tb/KidKTtJ °f the 
Unitarity Triangle. However, in contrast to the theoret- 
ically clean determination of the angle (3 in B° decays 
to charmonium final states 0, , the extraction of a in 
B° — > 7r + 7r - is complicated by the interference of tree 
and penguin amplitudes with different weak phases. The 
shift between a e g, from the measured B° — > ir + ir~ asym- 
metry, and a may be evaluated or constrained using mea- 
surements of the isospin-related decays B°(B°) — > 7r°7r° 
andP± ^t^tt 0. 

The CP-violating charge asymmetry for B^ modes, 
defined as 

_ \A\i-\A\* 

Acp= \ay+W () 

where A (A) is the B + (B~) decay amplitude, will de- 
viate from zero if the tree and penguin amplitudes each 
have different weak and strong phases. In the Standard 
Model the decay B^ — > 7r ± 7r° has only a tree ampli- 
tude contribution, so no charge asymmetry is expected. 
Both the rate and asymmetry of the decay P ± — » K^tt 
may constrain the value of the Unitarity Triangle an- 
gle 7. In particular, the ratio of B(B ± — ► K^tt ) and 
i3(B ± — > if°7r ± ) provides a lower bound for 7 [?]. The 
decay P ± — > K ± 7r° can also exhibit a significant charge 
asymmetry; different models for hadronic B decays pre- 
dict a range of values || . 

In this paper, we report on an observation of the decays 
P ± — > ^tt and B ± — » K ± tt° , a measurement of their 
CP- violating charge asymmetries, and a search for the 
decay B° -> 7r°7r°, using (87.9 ± 1.0) x 10 6 BB pairs 
collected with the BABAR detector. 

BABAR is a solenoidal detector optimized for the 
asymmetric-energy beams at PEP-II and is described 
in detail in Ref. |9J|. Charged particle (track) momenta 
are measured with a 5-layer double-sided silicon vertex 
tracker (SVT) and a 40-layer drift chamber (DCH) in- 
side a 1.5 T superconducting solenoidal magnet. Pho- 
ton (neutral cluster) positions and energies are measured 
with an electromagnetic calorimeter (EMC) consisting of 
6580 CsI(Tl) crystals. Tracks are identified as pions or 
kaons by the Cherenkov angle 6 C measured with a detec- 
tor of internally reflected Cherenkov light (DIRC). 

High efficiency for recording BB events in which one B 
decays with low multiplicity is achieved with a two level 
trigger with complementary tracking and calorimetry- 
based trigger decisions. BB events are selected using 
track and neutral cluster content and event topology. 



Candidate ir° mesons are reconstructed as pairs of pho- 
tons, spatially separated in the EMC, with an invariant 
mass within 3cr of the ir° mass. The resolution sigma is 
approximately 8 MeV/c 2 for high momentum ir°. Pho- 
ton candidates are required to be consistent with the ex- 
pected lateral shower shape, not be matched to a track, 
and have a minimum energy of 30 MeV. To reduce the 
background from false ir° candidates, the angle 9 1 be- 
tween the photon momentum vector in the tt° rest frame 
and the ir momentum vector in the laboratory frame is 
required to satisfy |cos# 7 | < 0.95. The 7r° candidates 
are fitted kinematically with their mass constrained to 
the nominal ir° mass. 

Candidate tracks are required to be within the track- 
ing fiducial volume, originate from the interaction point, 
consist of at least 12 DCH hits, and be associated with 
at least 6 Cherenkov photons in the DIRC. 

B meson candidates are reconstructed by combining 
a 7T° with a pion or kaon (h^) or by combining two 7r° 
mesons. Backgrounds arise from two sources: B — > pix 
decays in which one pion is emitted nearly at rest in 
the B frame so that the remaining decay products are 
kinematically consistent with a B ± — ► n^n or B° — > 
7r°7r° decay, and e + e~ — > qq{q = u,d,s,c) events where 
an or ir° from each quark randomly combine to mimic 
a B decay. 

Both backgrounds are separated from signal using 
the kinematic constraints of B mesons produced at the 
T(4S). The first kinematic variable is the beam-energy 
substituted mass toes = \J ( s /2 + P; ■ Pb) 2 /Ef — p^, 
where y/s is the total center-of-mass (CM) energy. 
(Ei,pi) is the four-momentum of the initial e + e~ sys- 
tem and pb is the B momentum both in the laboratory 
frame. The second variable is AE = Eb — \/s/2, where 
Eb is the B candidate energy in the CM frame. The 
pion mass is assigned to all /i ± candidates for the AE 
calculation. 

The .B — > p ± tt° background to B° — > 7r 7r° is reduced 
by only using candidates with \AE\ < 0.2 GeV. Re- 
maining B ± — > p^n background is further suppressed 
by removing candidates in which the additional ^ is 
identified. The track that gives a ir^n invariant mass 
and Toes of the ^ifiir combination most consistent 
with the p and B mass is selected. Requirements on 
the resulting 7r ± 7r° invariant mass and on the AE of 
the 7r ± 7r°7r° combination remove roughly 50% of the re- 
maining B ± — > p ± 7r° background, with 93% efficiency for 
B° -> 7r°7r°. Only (0.40 ± 0.04)% of B ± -> P ± tt° decays, 
and a negligible fraction of nonresonant B ± — > 7r ± 7r°7r° 
decays, remain after all cuts. For S ± — > k^ir the 
B — > pir background is suppressed by selecting candi- 
dates with —0.11 < AE < 0.15 GeV. 

The jet-like qq background is suppressed by requiring 
that the angle S between the sphericity axes of the B 
candidate and of the remaining tracks and neutral clus- 
ters in the event, in the CM frame, satisfy |cos# s | < 



5 



TABLE I: The results for both B — > h ± ir° and B° — > 7r°7r° are summarized. The number of B candidates N, total detection 
efficiencies e, fitted signal yields Ns, significances S, charge-averaged branching fractions B, asymmetries A, and 90% C.L. 
asymmetry limits are shown. Errors are statistical and systematic respectively, with the exception of e whose error is purely 
systematic. The upper limit for the B° — > 7r°7r° branching fraction corresponds to the 90% C.L., and the central value is shown 
in parentheses. 



Mode 


N 


e(%) 


N s 


S(a) 


S(10~ 6 ) 


A 


A(90%C.L.) 




21752 


26.1 ± 1.7 
28.0 ±2.0 


125±2? ± 10 
239±22 ± 6 


7.7 
17.4 


5.5±o% ± 0.6 
12.8ii? ± 1.0 


— o.o3±g:if ± 0.02 

-0.09 ±0.09 ±0.01 


[-0.32,0.27] 
[-0.24,0.06] 


7T%° 


3020 


16.5 ± 1.7 


90+IO +8 
zo -9 -4 


2.5 


< 3.6 (i.6±g;S ±S:i) 







0.8(0.7) for B ± -> h ± ir a (B° -► 7r°7r°). Also, we require 
toes > 5.2GeV/c 2 . The number of B ± -> h ± Tr° and 
candidates satisfying these requirements and 
the estimated efficiencies, obtained from simulated data, 
are shown in the first two columns of Table [I] The sim- 
ulation has been tuned to reproduce the observed track 
and 7T° efficiencies. The error in the estimated efficiency 
is dominated by the 5% systematic uncertainty in the 
single 7T° reconstruction efficiency. 

The number of signal B candidates is determined in an 
extended unbinned maximum likelihood fit. The proba- 
bility Vi (xj] di) for a signal or background hypothesis is 
the product of probability density functions (PDFs) for 
the variables Xj given the set of parameters di . The like- 
lihood function is given by a product over all N events 
and the M signal and background hypotheses: 



M 



N 



C = exp I - ^ nA Yl 



i=i 



' M 

E 

.i=l 



(2) 



For B^ — > /i ± 7r° the probability coefficients are Ni = 
|(1 — qjAi)rii, where qj is the charge of the track h and 
the fit parameters n, and A4 are the number of events 
and asymmetry for the four 7r + 7r° and K + ir° signal and 
background components. For B° — > 7r 7r° the coefficients 
are Ni = rij where the three rii are the number of signal 
candidates, B^ — ► p^n background and qq background. 
Monte Carlo simulations are used to verify that the like- 
lihood fits are unbiased. 

The variables Xj used for B ± — > ^tt are hies, AE, 
the Cherenkov angle 8 C of the /i ± track, and a Fisher 
discriminant T . The Fisher discriminant is given by an 
optimized linear combination of XaP* an( ^ cos#i| 2 
where pi is the momentum and 9i is the angle with re- 
spect to the thrust axis of the B candidate, both in the 
CM frame, for all tracks and neutral clusters not used to 
reconstruct the B meson. 

The PDFs for tubs, AE, 6 C , and T for the background 
are determined using data, while the PDFs for signal are 
found from a combination of simulated events and data. 
The jties distribution for background is modeled as a 
threshold function [lo| . whose shape parameter is a free 
parameter of the fit. The AE distribution for background 



is modeled as a quadratic function whose parameters are 
determined from the tties sideband in data. The tties 
and AE distributions for signal are modeled as Gaus- 
sian distributions with a low-side power-law tail whose 
parameters are found with simulated events. The AE 
resolution is approximately 42 MeV based on simulated 
events and is confirmed by evaluating the resolution in a 
sample of S ± — > D°p ± (p ± — * ^n ) events with an en- 
ergetic 7T°. To allow for EMC energy scale variations, 
the mean of the AE PDF is a free parameter of the 
fit. To account for the use of the pion mass hypothe- 
sis, the mean of AE is shifted for the K ± ir° PDFs. The 
T distribution is modeled as a bifurcated Gaussian and a 
double Gaussian for signal and background respectively, 
whose parameters are determined for signal from sim- 
ulation and for background from hies sidebands. The 
difference of the measured and expected values of 9 C for 
the pion or kaon hypothesis, divided by the uncertainty 
on 6 C , is modeled as a double Gaussian function. A 
control sample of kaon and pion tracks, from the decay 
D* + — > D°ir + , D° — * AT~7r + , is used to parameterize ag c 
as a function of the track polar angle. 

The variables Xj used for B° — > ir°ir° are mES) AE, 
and another Fisher discriminant J-q-. The Tt combines 
T with information from the B tagging algorithm de- 
scribed in Ref. The tagging algorithm uniquely clas- 
sifies events according to their lepton, kaon, and slow 
pion (from D* + — > -D°7r^ ow ) content, using all tracks in 
the event. Nine event classes, in decreasing order of their 
background rejection, contain the following: a high mo- 
mentum electron and a kaon, a high momentum muon 
and a kaon, a high momentum electron, a high momen- 
tum muon, a kaon and a slow pion, a well identified kaon, 
a slow pion, any kaon, or none of the above. These event 
classes are assigned an index, which is a new discriminat- 
ing variable, and is combined with T into a second Fisher 
discriminant Tt, optimized using simulated events. 

The toes distribution for qq background is param- 
eterized by the same threshold function used in the 
_B ± — * h^n analysis, where the shape parameter is de- 
termined from data with | cos(9 s | > 0.9. The AE distri- 
bution for qq background is modeled as a quadratic poly- 
nomial with parameters found from on-resonance data 
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B — > tt^tt (top) and B — > if tt (bottom), for candidates 
that satisfy optimized requirements on probability ratios for 
signal to background based on all variables except the one 
being plotted. The fraction of signal events included in the 
plots is 24% (mss) and 35% (AE) for t^tt , and 53% (m ES ) 
and 48% (AE) for K ± ty°. Solid curves represent projections 
of the complete maximum likelihood fit result; dotted curves 
represent the background contribution. For the B — > it ± tt 
AE distribution, the dotted curve shows the qq background 
and the small B — ► K ± tt° cross- feed; the dashed curve in- 
cludes the B — > pn background as well, so is the sum of all 
backgrounds. 



in the toes sidebands and off-resonance data. The Toes 
and AE variables in both B° — > tt°tt° and i? ± — > p^n® 
are correlated, so a two dimensional PDF derived from a 
smoothed simulated distribution is used. The AE res- 
olution is approximately 80 MeV. The distribu- 
tion for qq, B ± — > /9 ± 7r°, and B° — > tt°tt° is modeled 
as the sum of three Gaussians. For qq the parameters 
are found using both toes sideband and off-resonance 
data. For B° — » 7r°7r° and B ± — ► p ± ir Q the param- 
eters are found using a sample of fully reconstructed 
B" -> D^rar (n = 1, 2, 3) events. 

The decay B ± — > p*^ has not been observed; Ref. [ll| 



set an upper limit of B(B d 



p 7T ) < 4.3 x 10" 



at 



90% C.L. based on a measured central value of B(B^ — > 
p ± 7r°) = 2.4 x 10~ 5 . Therefore we fix the number of 
S ± — > p ± 7r° events in the fit to n p7r o = 8.4, based on this 
central value, and evaluate the systematic uncertainty of 
allowing n p ^o to vary from 4.2 to 15 events. 

The results of the maximum likelihood fits are summa- 
rized in Table |U Distributions of some of the variables 
used in the fits are shown in Figs.nand[21for B ± — > h^ir 
and B° — > 7r°7r°, respectively. The data shown are for 



events that have passed a probability ratio cut optimized 
to enhance the signal to background fraction. The like- 
lihood function for B° — > 7r°7r° is shown in Fig. 01. The 
statistical errors on the number of events are given by the 
change in signal yield n.j that corresponds to an increase 
in — 21n£ of one unit. The systematic uncertainty in the 
likelihood fit is estimated by varying the PDF parame- 
ters by their statistical errors or by comparing the result 
with an alternate parameterization. 

For B ± — > 7T ± 7r°, the dominant systematic uncertainty 
is due to the J- PDF for signal (±6.2 events) and back- 
ground (±7.6 events) PDFs, while for B ± — > K ± ir° it is 
due to the toes PDF for signal (j^'g events). System- 
atic uncertainties on the CP asymmetries are evaluated 
from PDF parameter variations and the upper limit on 
intrinsic charge bias in the detector (1.0%). 

For B° — > 7r°7r°, systematic uncertainties from the 



f+7.5 



PDFs are due to the Tr PDF for qq background (_ 2 a 
events), the toes PDF for qq background (~t{_i events), 



Ad- 



and the AE PDF for qq background (1" !2 events), 
ditional systematic uncertainties for B° 
from uncertainty in the EMC energy scale (iti 'f events), 
•° rejection cut (±1.3 events), and uncer- 



thc B ± -> p-TT" 
tainty in the assumed B — > p ±r K K} branching fraction 
(li'g events). The significance of the event yield, also 
listed in Table [I] is evaluated from the square root of the 
change in — 21n£ with the signal yield fixed to zero. The 
upper limit for B° — > 7r°7r° is evaluated by finding n n o^o 
where J™" " C(n)dn/ J °° C(n)dn = 0.9. For both sig- 
nificance and upper limits, systematic uncertainties are 
included with a worst case assumption for efficiencies and 
PDF variations. 

We observe B(B ± -► 7 r ± 7r ) = (5.5tJ;g ± 0.6) x 10" 6 , 
with a statistical significance of 7.7a from zero. This 
result is consistent with several prio r measurements re- 
porting evidence for this decay |l2l fl3l Q|. We mea- 
sure B(B ± -> ^tt ) = (12.8t\i ± TO) x 10" 6 . No 
evidence of direct CP violation is observed. Our limit 
B(B° — > 7r°7r°) < 3.6 x 10 -6 improves upon prior re- 
sults 0,^3- Removing correlated systematic uncertain- 
ties from luminosity and 7r° efficiency, we bound the ratio 
B(B° -> Tr°ir°)/B(B ± -> ^±77°) < 0.61 at a 90% confi- 
dence level. Assuming isospin relations for B — ► 7T7t 0, 
this corresponds to an upper limit of |a e ff — ot\ < 51°. 

We are grateful for the excellent luminosity and ma- 
chine conditions provided by our PEP-II colleagues, and 
for the substantial dedicated effort from the computing 
organizations that support BABAR. The collaborating 
institutions wish to thank SLAC for its support and 
kind hospitality. This work is supported by DOE and 
NSF (USA), NSERC (Canada), IHEP (China), CEA and 
CNRS-IN2P3 (France), BMBF and DFG (Germany), 
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FIG. 2: The a) m ES , b) AE, and c) Tr distributions for 
B° — + 7r°7r° are shown, for candidates that satisfy optimized 
requirements on probability ratios for signal to background 
based on all variables except the one being plotted. The frac- 
tion of signal events included in the plots is 20%, 20% and 
63% for niES, A_E and Tr , respectively. The dotted lines 
show the PDF projections for both qq and B — > p ± ty° back- 
ground, while the solid lines are the PDF projections for sig- 
nal plus background. The ratio — 21n (H/Cmam) is shown in 
d) where the dashed line is for statistical errors only and the 
solid line is for statistical and systematic errors, as applied 
for the calculation of significance. 
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